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ABSTRACT 

We consider constraints on the amplitude of mass fluctuations in the universe, erg, derived from 
two simple observations: the present number density of clusters and the amplitude of their correlation 
function. Allowing for the possibility that the primordial fluctuations are non-gaussian introduces a 
degeneracy in the value of as preferred by each of these constraints. However, when the constraints 
are taken together this degeneracy is broken, yielding a precise determination of cg and the degree of 
non-gaussianity for a given background cosmology. For a flat, f2m = 1 universe with a power spectrum 
parameterized by a CDM shape parameter F = 0.2, we find that the perturbations are consistent with a 
gaussian distribution with erg = 0.49tQ 07 (95% limits). For some popular choices of background model, 
including the favored low matter density models, the hypothesis that the primordial fluctuations are 
gaussian is ruled out with a high degree of confidence. 

Subject headings: Cosmology: observations and theory 
structure of the universe. 



galaxies: clusters: general — large-scale 



1. INTRODUCTION 

Most studies of structure formation in the universe as- 
sume that the primordial density perturbations are gaus- 
sian. Standard inflationary theories predict gaussian per- 
turbations, and the central limit theorem tells us that 
any theory involving the superposition of many random 
processes will give rise to approximately gaussian fluctu- 
ations. However, many well motivated theories predict 
non-gaussian initial conditions, including topological de- 
fect theories (Kibble 1976) and certain forms of inflation 
(Peebles 1983,1997). So far, no convincing observational 
evidence has been found to confirm or refute the gaussian 
hypothesis. 

Clusters of galaxies, however, provide us with a unique 
probe of possible non-gaussianity. Being the most mas- 
sive collapsed structures, they correspond to rare peaks 
in the primordial density field, so their statistics respond 
very sensitively to non-gaussianity in the initial matter 
distribution. In addition, an analysis of the formation 
of clusters from given initial conditions requires primar- 
ily gravitational physics, and is largely free of complica- 
tions from star formation and feedback. Finally, although 
the gravitational evolution which must be modeled is non- 
linear, well-tested analytical approximation schemes exist 
for studying the statistics of the resulting cluster distribu- 
tion. 

In this work we adapt these analytic schemes to the case 
of non-gaussian fluctuations, and use them to make pre- 
dictions for two simple observations, the number density 
of clusters and the amplitude of cluster correlations. We 
show that the possibility of non-gaussianity introduces a 
degeneracy in the value of CTg (where an is the rms over- 
density in a sphere of radius R /i~^Mpc) preferred by each 
of these observations. However, combining the two con- 
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straints breaks the degeneracy, allowing a precise deter- 
mination of (Tg and the degree of non-gaussianity in the 
universe. 

2. MODEL 

We describe the primordial density field by specifying 
both the power spectrum P{k) and the probability den- 
sity function (PDF) pr{5) of the fractional over-density S 
after a top-hat smoothing on scale R. We make the sim- 
plifying assumption that only the rms level of fluctuations 
changes as a function of scale, that is pn (6) = P{5/(jji) /gr 
where P{y) is a rescaled PDF, whose distribution has zero 
mean and an rms value of 1 . This assumption should hold 
reasonably well for topological defect theories, whose per- 
turbations are scale invariant over a range of scales, and 
more generally, will not affect our results provided the form 
of the PDF does not change significantly over the narrow 
range of scales relevant to cluster formation. 

The Press-Schechter (1974) formalism (hereafter PS) al- 
lows us to compute the number density of clusters within 
a given background cosmology for gaussian fiuctuations. 
The derivation suggests a simple way to adapt the PS for- 
mula for theories with a general PDF (as done by Chiu, 
Ostriker & Strauss 1997): The probability that a given re- 
gion of space lies inside a collapsed region with pre-coUapse 
radius larger than R is taken to be J^^ pii{6)dS, where Sc 
is a critical threshold for collapse. Differentiating with 
respect to R and dividing by the pre-coUapse volume we 
obtain the number density n{R) of collapsed objects with 
pre-collapse radius R, where 
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As in the gaussian case (Bond et al. 1991), we have mul- 



2 



tiplied by a correction factor / = 1/ Jp°° P{y)dy, to ensure 
that the entire mass of the universe is accounted for. We 
integrate the above expression to obtain cumulative num- 
ber counts A^>_R for all objects with a pre-coUapse radius 
greater than R: 

where ydR) = Sc/crn. We convert to cumulative number 
counts as a function of mass and temperature using the 
method of Eke, Cole & Frenk (1996). 

The peak-background split, first suggested by Kaiser 
(1984), provides a simple scheme for estimating the bias, 
and hence the correlation amplitude, of collapsed objects. 
Consider the manner in which the addition of a long wave- 
length background fluctuation 6b{x) modulates the local 
number density of collapsed objects. The threshold for 
critical collapse due to short wavelength fluctuations is 
modified locally from Sc to dc — '5b(x), so the number den- 
sity of collapsed objects is modified from its background 
vaiue N>R to 7V>fl(x) = (1 + Sbi-x.))NyR\s^^s^^s^{^). The 
factor 1 + (5f,(x) comes from the linear growth of perturba- 
tions. Taylor expanding to first order in db we obtain 

iV>flXx) = N^r{1 + 5b{^)) - Sb{^)^^. (3) 

dOc 

The large scale cluster-cluster correlation function is given 
by ^cc{x) = ((5c(x)5c(0)) where x = |x| and (5c(x) = 
(■^>fl(x) — Nyfi)/ NyR is the fractional overdensity in 
the cluster distribution. Defining bias b via ^cc(a^) — 
6^^(x) where ^(x) = {Sb{x)Sb{0)) is the correlation func- 
tion of the dark matter on large scales, we find b — 
l-{dNyR/d6c)/N>R. 

Substituting using equation |^ yields the bias: 

PMR)]/{R^aR)^3j^dR'PMR')]/iR'^aR,) 
dR'P[y,{R')]y,{R') {-daR,/dR') /{R'^gr,) ' 

(4) 

For each of our models we specify the background cos- 
mology in terms of and Oa- We then parameterize the 
power spectrum P{k) of matter fluctuations as a CDM 
spectrum (Bardeen et al. 1986, eq. G3) with shape pa- 
rameter r (where T ~ VLh) and normalization specified 
by CTg. Over the range of scales relevant to cluster for- 
mation this form gives a useful fit to power spectra in a 
variety of structure formation scenarios. We choose the 
PDF P{y) to be either a or a log-normal distribution, 
translated and renormalized to have mean zero and stan- 
dard deviation 1. In each case there is one free parameter 
which allows us to dial the amount of deviation from gaus- 
sianity, which we quantify in terms of a new parameter 
T = V^/°° P{y)dy/J^ e-y'/^dy, where T is the prob- 
ability of obtaining a three or higher standard deviation 
peak for the PDF in question relative to that for a gaus- 
sian PDF. To obtain log-normal and distributions with 
T less than 1, we reflect the PDF by making the transfor- 
mation P{y) P(— y). Some sample PDFs, along with 
the value of T, are illustrated in Fig. 1. Our models are 
thus specified by five parameters: flm, ^A, F, as and T. 

The PS formula is based on analytic results for the col- 
lapse of spherically symmetric perturbations. In the case 



of gaussian fluctuations, it provides remarkably good flts 
to the multiplicity function of cluster scale objects (Gross 
et al. 1997). For non-gaussian fluctuations, one may worry 
that the spherical collapse picture is less applicable, partic- 
ularly in the case T <C 1, where the highest density objects 
will tend to form from matter swept out of low density re- 
gions. As for cluster correlations, our result for the cluster 
bias in the gaussian case is identical to that of Mo & White 
(1996), which has been shown by Mo, Jing & White (1996) 
to give a good fit to the amplitude of cluster correlations 
in N-body simulations. For highly non-gaussian fluctua- 
tions, problems could arise in the separation of the long 
and short wavelength fluctuations, which are not necessar- 
ily uncorrelated. However, if most clusters form from rare 
peaks in the initial density field, and if the long and short 
wavelength modes are nearly uncorrelated, then there is 
no reason to expect the modified PS or peak-background 
split formalisms to lead to errors. Our predictions seem to 
be born out qualitatively by the results of published sim- 
ulations of non-linear gravitational clustering from non- 
gaussian initial conditions (Park, Spergel & Turok 1991; 
Weinberg & Cole 1992; Borgani et al. 1994). However, a 
detailed quantitative comparison with numerical simula- 
tions remains to be carried out. 

3. RESULTS 

We compare our predictions with two sets of observa- 
tional data. Firstly, we compare predictions for the clus- 
ter temperature function at z = 0.05 with observations 
computed from the catalogue of nearby clusters compiled 
by Henry & Arnaud (1991). We compute the cumula- 
tive number counts at two values of the temperature, and 
estimate the covariance matrix by bootstrap resampling. 
For each resampling, we add gaussian random errors to 
the temperature consistent with the quoted observational 
uncertainty, and an additional 10% systematic error to al- 
low for uncertainties in the M vs. T normalization (Eke, 
Navarro & Frenk 1998). The observational data, along 
with 2(T-errorbars for our two chosen temperatures, are 
shown in the top panel of Fig. 2. 

Secondly, we compare with the amplitude of the cluster- 
cluster correlation function, quantified by ro, the value of 
r for which ^ccif) = 1- Croft et al. (1997) (hereafter 
C97) have computed tq as a function of mean cluster sep- 
aration d for subsamplcs of the APM cluster survey with 
different lower bounds on the cluster richness. Their re- 
sults are shown in the bottom panel of Fig. 2. We com- 
pare model predictions with just the d — 30/i~^Mpc data 
point, since this point is derived from the largest number 
of clusters, and the points are highly correlated. We use 
ro = 13.2 ± 0.91 ± 1.0/i~^Mpc, where we have corrected 
the correlation amplitude for redshift space distortions by 
subtracting l/i~^Mpc, as suggested by the simulations in 
C97. The first errors are their 2a statistical errors and the 
second are 2a errors we have added to account for system- 
atic uncertainties in the redshift space correction. For each 
model we use the techniques described above to compute 
the cluster correlation length ro at z = for clusters with 
pre-collapse radius determined by their mean separation 
d = N^l/^ = 30/i-iMpc. 

In Fig. 2 we compare observations and predictions for 
a set of four models, each with flm = 1, ^^a = and 
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r = 0.2. Models 1, 2 and 3 have a normalization given by 
(Tg = 0.5. Model 1 has gaussian fluctuations, model 2 has 
T = 0.0348 (a reflected distribution with 50 degrees of 
freedom), and model 3 has T = 16.3 (a with 1 degree 
of freedom). Model 4 has the same PDF as model 3, but 
a lower normalization, erg = 0.2. 

Prom the top panel of Fig. 2 we see that increasing T 
increases the number density of clusters of a given tem- 
perature. The reason for this is clear from Pig. 1: for 
the value of erg considered here, only peaks of height 
roughly 3tT and higher (shown by the heavily shaded re- 
gion under each curve) will collapse, since the ratio of 
the critical overdensity to the standard deviation o-g in 
a sphere which will collapse to form a typical mass cluster 
is Uc = Sc/cTg. = 1.69/0.5 ~ 3. So in this "standard" case, 
the parameter T directly counts the expected number of 
clusters. 

Now we consider the effect of changes in the PDF on 
the amplitude of the cluster correlation function. Prom 
the bottom panel of Fig. 2 wc sec that increasing T de- 
creases the correlation amplitude, for reasons which are 
also illustrated in Pig. 1. Adding a small background fluc- 
tuation modifies the local value of the critical overdensity 
for collapse, leading to a spatial modulation in the num- 
ber density of collapsed objects. The lightly shaded re- 
gions under the curves in Pig. 1 show the effect for each 
PDF of modifying from 3.0 to 2.5. Por the T > 1 
PDF which has a long flat tail, this modification produces 
only a small enhancement in the local number of objects 
capable of collapse, and hence there is only a small en- 
hancement in the clustering amplitude of those collapsed 
objects. Por the T < 1 PDP where the positive tail falls 
off very sharply, this modification greatly enhances the lo- 
cal number of collapsed objects, and consequently there 
is large enhancement in the clustering amplitude of those 
objects. 

Finally, we consider the effect on each of the above re- 
siilts of changing erg for a fixed PDF. Model 4 illustrates 
the effect of decreasing ag for the case T > 1. As ag is 
decreased, the number density of clusters decreases, the 
reason being that the number of peaks with high enough 
over-density to collapse is decreased. The top panel of 
Fig. 2 shows that for T = 16.3, cg = 0.2 provides a good fit 
to the temperature function data. In the bottom panel of 
Pig. 2, we see that as ag is decreased, the amplitude of the 
cluster-cluster correlation function also decreases. This is 
due to the fact that decreasing erg decreases the ampli- 
tude of the underlying dark matter correlation function, 
but leaves the bias roughly unchanged. Model 4 shows 
that for T = 16.3, erg = 0.2 provides a worse fit to the 
correlation function data than erg = 0.5. 

Thus we see that the two datasets we have chosen are 
complementary: variations in the parameter T introduce 
a degeneracy in the preferred value of erg for each dataset, 
but the degeneracies are "orthogonal". For the tempera- 
ture function, increasing T requires lower erg, while for the 
correlation amplitude, increasing T requires higher erg. So 
taken in combination, these data allow a precise deter- 
mination of both (Tg and T. We quantify this by com- 
puting values of the x^ statistic for each model relative 
to the two datapoints in the top panel of Fig. 2 and the 
d = 30/i~^Mpc data point in the bottom panel. In Fig. 3 



we show goodness of fit contours in the erg vs. T plane 
for log-normal PDF models with a range of background 
cosmologies and values of the parameter P. (We use three 
values of P which span the range allowed by measurements 
of the APM cluster power spectrum by Tadros, Efstathiou 
& Dalton 1997.) 

The dotted blue line (T = 1) shows the location in the 
T vs. erg plane of a gaussian PDP. The constraints on erg 
from present day cluster number abundances have been 
extensively investigated in this restricted case (White, Ef- 
stathiou & Frenk 1993; Eke et al. 1996; Viana & Liddle 
1996) and the results obtained here are in close agree- 
ment. In addition. Mo et al. (1996) have investigated the 
combined constraints of cluster correlation amplitude and 
number abundance measurements for the case of gaussian 
fluctuations. Our results along this "line of gaussianity" 
can be considered an update of their analysis, as we use 
more recent cluster correlation data (C97) with smaller 
observational errors. 

The results depend both on the background cosmology 
and the power spectrum shape parameter P. For low val- 
ues of flm, the value of erg preferred by the temperature 
function data increases. This is just as observed in the 
gaussian case, and is due to the fact that for low values 
of rim, clusters of typical mass form from spheres with 
a larger pre-coUapse radius R, for which an is less than 
erg. The value of ag preferred by the correlation amplitude 
data, however, is almost independent of Qm and Q\. This 
is because the correlation amplitude is measured for clus- 
ters with a fixed value of d, and consequently for clusters 
with the same pre-coUapse radius (not the same mass). 
Finally, as P increases, the value of erg preferred by the 
correlation amplitude also increases. This is because the 
corresponding change in the shape of the power spectrum 
decreases the amount of power on scales of order the clus- 
ter correlation length for a given value of ag. 

Results for the best fit values of erg and T from the 
conibin(xl datasets for the family of log-normal PDFs in 
the various background cosmologies are summarized in ta- 
ble 1. We have performed the same calculation for x^- 
parameterized PDFs and find very similar results. Thus 
we expect that we can infer results for other models pro- 
vided we work out the appropriate value of the parameter 
T (for instance T ~ 14 for textures - Park et al. 1991). In 
future work we will expand on the analysis presented here 
by including observations of the evolution of cluster abun- 
dance and making use of updated datasets at low redshift. 

4. CONCLUSIONS 

We have modified the Prcss-Schcchtcr and the peak- 
background split formalisms to compute the cluster num- 
ber density and cluster correlation amplitude in non- 
gaussian models. The best-fit value of erg and level of non- 
gaussianity depend on the choice of background cosmology, 
but wc have demonstrated how these two quantities can 
be constrained simultaneously. For a flat, f^A = 0, P = 0.2 
universe we find that the fluctuations are consistent with 
gaussian, with erg = 0.49lo;557 (95% limits). For low-f7„ 
models with F = 0.2, gaussian fluctuations are ruled out 
with a high degree of confidence, while non-gaussian fluc- 
tuations (T > 1) are strongly preferred. Taken in con- 
junction with accurate knowledge of the background pa- 
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rameters (which we are rapidly gaining from supernovae 
and CMB observations), this method can provide a power- 
ful constraint on the nature of density fluctuations in the 
universe. 

Some concerns remain regarding the validity of our mod- 
ification of the PS and peak-background split formalisms 
to the case of non-gaussian fluctuations, and some system- 
atic errors may not yet be accounted for. This issue should 
be settled by an analysis of n-body simulations with non- 
gaussian initial conditions. However, our formulae have 
been well tested for the gaussian case, so the result that 



gaussianity is ruled out for certain choices of background 
model is not subject to these uncertainties. At the very 
least, we have demonstrated the power of two simple and 
well measured datasets, the cluster temperature function 
and the amplitude of cluster correlations, to simultane- 
ously constrain the amplitude of mass fluctuations and 
the level of non-gaussianity in the universe. 
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supported in part by a grant from the NSF, and E.G. ac- 
knowledges the support of an NSF Graduate Fellowship. 
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Fig. 1. — Three PDFs. The heavily shaded region shows the area contributing to peaks of height 3a and above, while 
the lightly shaded region shows the additional contribution to peaks of height 2.5a and above. 
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Table 1 

Best fit values of as and T, with 95% confidence limits. 
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Fig. 2. — The cluster temperature function (black jagged line, top panel), and the cluster correlation length ro (dat- 
apoints, lower panel), with predictions for four models (f2„i = 1-0, f^A — 0.0, and F = 0.2, erg and T labeled in lower 
panel). 
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Fig. 3. — Confidence intervals (68% — dark band, 95% — light band) in the T vs. cts plane for the cluster temperature 
function (green band, top left to bottom right) and the cluster correlation length (red band, bottom left to top right). 
The blue dotted line corresponds to gaussianity [T = 1), and the panels are labeled with the values of f^m, ^^A, and F. 



